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Abstract 

The statistical properties of the multiphers of the absolute returns are investigated 
using one-minute high-frequency data of financial time series. The multiplier dis- 
tribution is found to be independent of the box size s when s is larger than some 
crossover scale, providing direct evidence of the existence of scale invariance in fi- 
nancial data. The multipliers with base a = 2 are well approximated by a normal 
distribution and the most probable multiplier scales as a power law with respect 
to the base a. We unravel that the volatility multipliers possess multifractal nature 
which is independent of construction of the multipliers, that is, the values of s and 
a. 
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1 Introduction 



It has been a long history that physicists show interests on financial mar- 
kets, which can be at least traced back to 1900 when Bachelier modeled stock 
prices with Brownian motions [1]. In the middle of last century, Mandelbrot 
proposed to characterize the tail distributions of income and cotton price fluc- 
tuations with the Pareto-Levy law and applied R/S analysis to investigate the 
temporal correlations in the evolution of stock prices [2] . In recent years since 
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the seminal work of Mantegna and Stanley [3], Econophysics has attracted 
extensive interest in the physics community. 

As an analogue to turbulence, many time series observed in the financial mar- 
kets are reported to possess multifractal properties [US], such as the foreign ex- 
change rate [HISllSllTllSll^lUIIll] , gold price [S] , commodity price [12] , stock price 
plITMllTSllMTTj . stock market index [TMM20]|2T]|22ll23]IMll2^^ to 
list a few. Extensive methods have been adopted to extract the empirical mul- 
tifractal properties in financial data sets, for instance, the wavelet transform 
module maxima (WTMM) [2^3U|31j and the multifractal detrended fiuctua- 
tion analysis (MF-DFA) [52] • A time series of the price fiuctuations possessing 
multifractal nature usually has either fat tails in the distribution or long- 
range temporal correlation or both [32] . However, possessing long memory is 
not sufficient for the precence of multifractality and one has to have a nonlin- 
ear process with long-memory in order to have multifractality [33]. In many 
cases, the null hypothesis that the reported multifractal nature is stemmed 
from the large fiuctuations of prices can not be rejected |34j . 

Here, we propose to investigate the multifractal nature of absolute returns 
of stocks based on the multiplier method, again, borrowed from turbulence 
[35|36|37ll38] . Our goal is to provide direct evidence of scale invariance in the 
distribution of the multipliers. The concept of multiplier was originally in- 
troduced by Novikov to describe the intermittency and scale self-similarity 
in turbulent fiows |3M|. The scale-invariant multiplier distribution is argued 
to be more basic than the standard /(«) curve [35)136] . In addition, it allows 
us to extract both positive and negative parts of the /(a) function with ex- 
ponentially less computational time and is more accurate than conventional 
box-counting methods [3SI3S]- 



2 Description of the data set 

We adopt a nice high-frequency data set recording the S&P 500 index to 
ensure better statistics in our analysis. The record contains quoted prices I{t) 
of the index, covering eighteen years from Jan. 1, 1982 to Dec. 31, 1999. The 
sampling interval is one minute. As usual, the nontrading time periods are 
treated as frozen such that we count merely the time during trading hours 
and remove closing hours, weekends, and holidays from the data. The size of 
the data set is about 1.7 million. 

The return r{t) over a time scale At is defined as follows 

r(t)=ln[/(t)]-ln[/(t-At)] , (1) 
whose absolute is a measure of the volatility. In this letter, the time scale is 
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At = 1 mill. We can construct an additive measure in the time interval [ti, ^2], 
which is the sum of absolute returns: 

K[tiM) = f:\r{t)\. (2) 

t=ti 

The quantity /x([ti,t2]) is actually a measure of the volatility on the time 
interval [^1,^2] [lO]- The time series \r{t) \ is partitioned into boxes of identical 
size s. Each of these mother boxes is further divided into a daughter boxes. 
The multiplier m is determined by the ratio of the measure on a daughter box 
to that on her mother box Therefore, the multiplier is dependent of a 
and s and can be denoted as when necessary. 



3 Scale invariant distribution 

Figure [T](a) presents the probability densities pa,s{'m) of the multiplier m for 
four different box sizes s = 30, 60, 120, and 180 with the same base a = 2. All 
curves are symmetric with respect to m = 0.5 such that P2,s('^) = P2,s(l — fn) 
by definition and close to Gaussian. The solid lines are the best fits to normal 
distributions, whose fitted standard deviations are a = 0.089 for s = 30, 
a = 0.073 for s = 60, a = 0.068 for s = 120, and a = 0.069 for s = 180, 
respectively. The corresponding r.m.s. of the fit residuals are 0.050, 0.068, 
0.050, and 0.068. Note that the mean /i = 1/2 is fixed in the fitting procedure. 
It is evident that a decreases in regard to s and tends to a constant for large s. 
This phenomenon is further manifested by Fig. [U^b), which plots the sample 
standard deviation a of the multipliers as a function of the box size s for 
different base a. The inset shows the loglog plots of a against s. 




Fig. 1. (color online), (a) Empirical probability density functions of the volatility 
multipliers with the base a = 2 for different box sizes. The solid lines are fitted nor- 
mal distributions, (b) Dependence of the sample standard deviation ex as a function 
of s for different bases. The inset shows the loglog plots of a against s. 
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One can see that there are two regimes in the a versus s relation: a decays 
as a power law for small s and saturates to a constant for large s. Roughly- 
speaking, the crossover values of s are the following: Sx ~ 80 for a = 2, 
Sx ^ 100 for a = 3, Sx ~ 120 for a = 4, Sx ~ 140 for a = 5, Sx ~ 150 
for a = 6, and Sx ~ 150 for a = 7, respectively. In other words, the sample 
variance o"^ in the case of a = 2 has the fastest convergence rate to a constant. 

Figure [2t^a) shows the empirical probability density functions pa^s{fn) of the 
multipliers for different bases a = 2, 3, and 5 and different box sizes s = 150, 
210, and 300. For each CL, Pa,sij^^ remains invariant in respect to s when s > Sx. 
In other words, there is a scaling range in which the volatility multiplier is 
scale invariant, whose distribution is independent of s. We shuffled the return 
series and found that the multiplier distributions are not scale invariant and 
the scaling range of s disappears, as illustrated in Fig.[2](b). Therefore, pa,s{jn) 
can be reduced to Pai^n) in the scaling range. The shuffling test shows that 
long memory in the volatility plays an essential role in the appearance of scale 
invariance. 
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Fig. 2. (color online), (a) Scale invariant distributions of volatility multiplier rria^s 
for different bases a = 2, 3, and 5 and different box sizes m = 150, 210, and 300. (b) 
Scale invariant distributions of volatility multiplier m^^s for the shuffled time series. 

The most probable multipliers are also investigated in this work. For a given 
s, the most probable multiplier mmax with base a is estimated such that 
Pa,s ('^max) = ^^^[Pa,s{^)]- Figure [3t^a) presents the loglog plots of mmax ver- 
sus a for different values of s. One can observe that the data points for different 
s collapse on a single line, showing a power-law dependence 

1^ma,x ^ (1 ^ , (3) 

where (3 = 1.10 ± 0.02 for s = 150, (3 = 1.10 ± 0.01 for s = 210, and (3 = 
1.09 ±0.01 for s = 300, respectively. Intuitively, since a mother box is divided 
into a daughter boxes, the sum of the a multipliers is one and the multipliers is 
expected to aggregate around 1/a. However, it is noteworthy that this power- 
law dependence is nontrivial, which does not hold in turbulence [36)1411142] . In 
Fig. Mh) is shown Pa,s ('^max) as a function of s for different a. It is evident 
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that p(mmax) 

increases with s and then approaches a plateau when s > Sx- 
Figure [3] further verifies the scale- invariant nature of the volatility multiplier. 



10 



10" 



S 



10" 



10 



10 



10 





s = 150 




5 = 210 




s = 300 


(a) 





10 



14 




a 



Fig. 3. (color online), (a) Power law dependence of mmax versus a for s = 150, 210, 
and 300. (b) The saturation behavior of p(minax) iu regards to s for different a. 



4 Multifractal analysis 



It is shown that, for any two bases a and b, the density functions are related 
through Mellin transform in the following form [36] 

[MK(m„)}]Vi- = [Mip.im,)}]'/'^" , (4) 

where M stands for Mellin transform. Equivalently, we have 

In /o mlpa{ma)dma _ In Jq mlph{mb) drrib 



In a In 6 



(5) 



The scaling exponent r(g) of the moment of can be obtained as follows 
13511361 



r(g) = -Do - , (6) 

ma 

where Dq is the fractal dimension of the support of the measure. In the cur- 
rent case, we have Dq = 1. Note that the use of the Mellin transform may 
indeed appear natural in the framework of scaling and power-like functions, 
for instance, in the analysis of Weierstrass-type functions [13] • 

The local singularity exponent a and its spectrum /(a) are related to T{q) 
through Legendre transforms: a{q) = T'{q) and /(a) = qa{q) — r^q). It follows 

that BSESHIHSHg ^^^^ ^ ^^^^^ 

{rria) ma 

and 

(mg)ln(mg) - {mlXnml) 
^^"^ = • 
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Equations (jMH]) predict that each of these characteristic multifractal argu- 
ments for fixed q converges to constant in the scahng range and are indepen- 
dent of the base a as well. 



In order to test this prediction, we have to calculate first the scaling func- 
tion r(g) which requires that the integrand m1pa{ma) converges for a given 
q imiH]. We have investigated mlpa{ma) for different values of q, s, and a. 
A typical dependence of mlpa{ma) as a function of is shown in Fig. IHfor 
different values of q with fixed box size s = 210 and base a = 2. The integrand 
diverges for large rria when q is larger than 6. We shall nevertheless investigate 
scaling functions for g ^ 8 for comparison. Moreover, Fig. H] indicates that the 
integrand diverges when g ^ — 1 and the associated negative moments do not 
exist. This is a direct consequence of the fact that Pa{0) ^ 0. Indeed, there 
are time moments when the local returns are zero so that the probability den- 
sity at = is apart from zero, i.e., Pa(0) 7^ 0. Approximately, for a small 
number 5, p{m) = p(0) is a constant for m < 6. Posing // m'^p{m)dm = C, 
we have mflp{m)dm = Jq m'^p{m)dm + C = p{0)-^m'^~^^ \ q + C. Therefore, 
we have g > — 1, which is quite analogous to the situation in turbulence 
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Fig. 4. Numerical integrand m'^p{m) as a function of m for different values of q with 
s = 210 and a = 2. 

We now present in Fig. [5] the results of the multifractal analysis for q varying 
from —1 to 8. In Fig.[5t^a) is shown the dependence of the scaling function r(g) 
upon the box size s for different values of q and a. It is evident that, for large 
s, T{q) is independent of s for every q under investigation. The r(g) function 
reaches constant faster for small a. These results are in excellent agreement 
with the theoretical predictions. Figure El^b) shows the dependence of the 
singularity spectrum f{a) on the box size s for different values of q and a. 
Again, we witness a range of scale invariance in which /(a) is independent of 
s. What needs to be emphasized is that, for q = 8, the three /(a) curves with 
different a do not converge due to bad statistics as shown by the left-middle 
panel of Fig. HI We plot the three scaling functions r(g) for a = 2, 3, and 5 
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with respect to q in Fig. [5]^c). The error bars are estimated as the standard 
deviation over different s. Except for large g, the three r(g) curves collapse 
on a single nonlinear curve. In addition, Fig. E](d) shows the three singularity 
spectra /(a) in respect to the local singularity exponent a for the three bases. 
Again, the three curves collapse remarkably on a single curve when q is not 
too large. Both Fig. [5](c) and Fig. strongly indicate that the volatility 
multiplier possesses multifractal nature. 
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Fig. 5. (color online). Multifractal analysis of absolute returns, (a) Dependence of 
r(g) with respect to s for different values of q and different bases a = 2, 3, and 5. 
(b) Dependence of /(a) with respect to s for different values of q and different bases 
a = 2, 3, and 5. (c) Scaling function r(g) obtained from Eq. ([6]) for different bases 
a = 2, 3, and 5. (d) Singularity spectrum f{a) obtained from Eq. ([7]) and Eq. ([8|) 
for different bases a = 2, 3, and 5. 

An important feature of multifractals is the possible existence of negative di- 
mensions in the multifractal spectrum, that is, /(a) < for large or small 
a [50]. Negative dimensions are more common if the multiplier distribution 
is continuous [36]l44f45l|l6] . Figure [5](d) also shows that there are negative di- 
mensions for large q especially when a is large. However, we should be cautious 
that the part of /(a) < might be an artifact of bad statistics for large q 
and a, as shown by the right panel of Fig. H] and Fig. [5](b) as well. Since the 
multifractal functions are more reliable statistically for a = 2, we argue that 
there is no negative dimension for g ^ 8. More data are required to investigate 
higher order moments and the issue of the existence of negative dimensions is 
still open. 

In the development of Econophysics, the literature has witnessed increas- 
ing analogues between turbulence and finance. Multiplier analysis is a well- 
established method in the description of conservative quantities in turbulence 
[39|35f36ll37f38] . Novikov predicted that the multiplier distribution p{m) is in- 
dependent of the scale s as long as s is well inside the inertial subrange [39] . 
Our finding provides further evidence of the analogue between turbulence and 
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finance, except that the existence of negative dimensions that was reported in 
turbulence is not confirmed in our financial data. Moreover, the energy dissi- 
pation multiplier with a = 2 follows approximately a triangular distribution, 
which is much flatter than the normal distribution of volatility multiplier in 
this work. 



5 Conclusion 

In summary, we have employed the multiplier method to investigate the volatil- 
ity of high-frequency data of the S&P 500 index. The distribution of volatility 
multiplier is found to be independent of the time scale s for different a when 
s is larger than some crossover scale Sx- We unraveled that the volatility mul- 
tipliers are scale invariant and have multifractal nature, which is independent 
of the construction of the multipliers (characterized by s and a) in the scaling 
range. 
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